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The Living With a Star Ion Flux Model (IFM) is a radiation environment risk mitigation 
tool that provides magnetospheric ion flux values for varying geomagnetic disturbance levels 
in the geospace environment. IFM incorporates flux observations from the Polar and 
Geotail spacecraft in a single statistical flux model. IFM is an engineering environment 
model which predicts the proton flux not only in the magnetosphere, but also in the solar 
wind and magnetosheath phenomenological regions. This paper describes the ion flux 
databases that allows for IFM output to be correlated with the geomagnetic activity level, as 
represented by the Kp index. 


I. Introduction 

T7 ven outside of the radiation belts, spacecraft can be exposed to a significant proton flux environment. The 
-^—'source of the proton flux can be from the solar wind, magnetosheath, or magnetosphere regions of geospace. 
Episodic injections of plasma from the magnetotail during substorms and major magnetic storms can increase proton 
flux in the energy band of concern by orders of magnitude in the outer magnetosphere. Although at flux levels that 
are typically lower than seen in the magnetosphere, protons are also found outside of the magnetosphere in the dusk 
and dayside magnetosheath, or even upstream of the bow shock, since “leakage” across the magnetopause is one of 
the loss mechanisms for magnetospheric plasma. Energetic solar protons are another source of potentially damaging 
particles. Solar event protons are a concern not only while a spacecraft is in the solar wind, but they also pose a risk 
while the spacecraft is in the magnetosheath or in the outer magnetosphere. Energetic solar protons easily traverse 
the bow shock and magnetosheath with little variation in flux and can even penetrate the lower magnetic field 
regions of the outer magnetosphere. 
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Figure 1. Sample CRM output. Kp dependent proton flux (in units of #/cm2-s-sr-MeV) given by the model for 
protons between 100 keV and 200 keV (adapted from Blackwell et al., 2003). 
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IFM addresses the need for a low energy flux model of the Earth’s magnetosphere that can be used in spacecraft 
design and mission analysis. This is an empirical engineering ion flux model applicable over a range of geocentric 
radial distances of approximately 2 Earth radii (Re) to the magnetopause within a maximum distance of 30 Re in the 
direction of the magnetotail. The model provides differential energy flux over a range from 58.1 keV to 1560.8 keV 
in 8 logarithmic energy steps. A key feature of the model is ion flux values traceable to satellite observations of the 
ion environment. Data is binned in spatial volume elements and sorted to yield Kp dependent flux estimates that are 
reported for percentile values (e.g., 50 %, 95 %) in the region. 

IFM is an alternative to standard environment models such as AP-8 which are not applicable for radial distances 
that include the outer magnetosphere where there are strong azimuthal asymmetries in charged particle flux. In 
addition, the modeling approach used in creating IFM includes not only the magnetosphere, but also the 
magnetosheath and near-Earth solar wind particle flux in the database, since these environments are available in the 
datasets used to derive the model. Thus, IFM provides an integrated charged particle environment definition model 
useful for a variety of applications including spacecraft design analysis, near real time environment modeling, and 
dynamic flux estimates for use in on-orbit operations support within 30 Re of the Earth. 

II. Technical Approach and Methodology 



Figure 2. Low energy (58.1 - 77.3 keV) proton flux as a function of the Kp magnetic activity index. 
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Figure 3. Medium energy (227.5 - 341.6 keV) proton flux as a function of the Kp magnetic activity i 


The development of IFM has heavily leveraged experience gained in creating the ion flux model for the 
Chandra X-ray observatory. The Chandra Radiation Model (CRM) was developed in response to the need for a tool 
to predict the 100-200 keV proton flux along its orbit [ Blackwell et al., 2000, 2003]. NASA standard trapped proton 
AP-8 [Sawyer and Vette, 1976] and electron AE-8 [Teague and Vette, 1974; Vette, 1991] models are used by the 
Chandra program to determine the mean locations of the very energetic radiation belts but are not adequate to 
address the issue of low energy protons in the outer radiation belts. CRM is the first engineering-level ion 
environment model for the outer magnetosphere. CRM is designed for use both as a scheduling tool for planning 
science observations for periods up to three weeks and for a real-time environment model for estimating low energy 
proton environments. Proton flux values as a function of Kp geomagnetic activity index from the CRM Version 2 
are shown in Fig. l.IFM is an empirical ion flux model of the magnetosphere using satellite 
observations of ion flux. Input to the model is location in space, date and time of year, and appropriate 
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Figure 4. High Energy (1560.8 - 3005.4 keV) proton flux as a function of the Kp magnetic activity index. 

interplanetary environment conditions (e.g., interplanetary magnetic field, solar wind density and velocity) and 
geomagnetic indices (e.g., Dst, Kp). The software returns values of the ion flux for user selected percentile levels 
(e.g., the maximum flux value that would be predicted to occur 50% or 95% of the time). 


IFM incorporates code from two separate space environment models to calculate magnetopause and bow shock 
boundary locations. The magnetopause model is taken from the Tsyganenko geomagnetic field model [ Tsyganenko , 
1995]. The bow shock model used is one adapted from a model by Bennett et. al. [Bennett, 1997]. The original 
version of CRMFLX provided ion flux values as a function of Kp but was limited both in its range of spatial 
application (-8 Re < Z G sm < +15 Re) and particle energy (100 keV to 200 keV protons) to meet specific 
requirements of the Chandra Program. 


2.1 Model Implementation 

Ion flux environments in the outer magnetosphere are complex and variable. Traditional techniques used to 
create trapped particle models (e.g., simple B-L flux mapping) are not applicable in magnetospheric regions where 
the geomagnetic field is highly perturbed and is significantly different than the dipole configuration. CRMFLX 
consists of a combination of analytical and database driven models, driven by proxy parameters of geomagnetic 
activity levels, K p and D st , that provide for correlation of magnetospheric particle flux with geomagnetic 
disturbances. 




Table 1. Geotail EPIC/ICS Energy Bands 


Channel 

Species 

Energy Band 

Sector 

Time Res. a 

Time Res. Database 13 

Name 


(keV/e) 

(deg) 

(sec) 

(sec) 

P2 

p+ 

58.1 

- 77.3 

22.5 

6 

288 

P3 

p+ 

77.3 

- 107.4 

22.5 

48 

288 

P4 

p+ 

107.4 

- 154.3 

22.5 

48 

288 

P5 

p+ 

154.3 

- 227.5 

22.5 

48 

288 

P4 

p+ 

107.4 

- 154.3 

22.5 

48 

288 

P5 

p+ 

154.3 

- 227.5 

22.5 

48 

288 

P6 

p+ 

227.5 

- 341.6 

22.5 

48 

288 

P7 

p+ 

341.6 

- 522.5 

22.5 

48 

288 

P8 

p+ 

522.5 

- 813.5 

22.5 

48 

288 

P9 

p+ 

813.5 

- 1560.8 

22.5 

96 

288 

P10 

p+ 

560.8 

- 3005.4 

22.5 

96 

288 


a Time resolution of original data. 

b Time resolution of spin averaged data obtained from Principle Investigator. 
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Figures 2-4 shows that there is good correlation between Kp and the ion flux in the outer magnetosphere. 
Unfortunately, the spatial regions that are sampled by spacecraft is very sparse during periods when the magnetic 
activity level is high (with correspondingly high ion flux). It is not possible then to have a model based strictly upon 
spacecraft measurement of the ion flux levels, while still having good coverage of all spatial regions of interest. 
Instead, the approach is a combination of database and analytical techniques that fill in the spatial gaps while at the 
same time maintaining a direct link to the satellite measurements. 

Ion flux measurements from the Geotail and Polar spacecraft were used to create the IFM database. Geotail is 
in a highly elliptical orbit that lies in the equatorial plane with a perigee of 8 Earth radii (Re) and apogee of 30 Re. 
This orbit makes Geotail particularly useful in creating an ion flux model for the outer magnetosphere, since it can 
directly sample the high flux plasma sheet region. While the equatorial orbit does allow for mapping along 
magnetic field lines to populate a larger region of the magnetosphere than what is directly measured, the Geotail 
data cannot provide information to populate the higher magnetic latitudes. Fortunately, the other principal data 


Table 2. CEPPAD/EPS Energy Bands 


Energy Species Energy Threshold Set 1 Energy Threshold Set 2 

Channel Min (keV) Mid(keV) Min (keV) Mid(keV) 


0 

H+ 

16.8 

18.9 

13.9 

15.6 

1 

HT 

21.2 

24.4 

17.5 

19.9 

2 

IT 

27.9 

32.4 

22.6 

26.2 

3 

IT 

37.5 

43.1 

30.3 

35.4 

4 

Hr 

49.6 

57.2 

41.4 

48.1 

5 

FT 

65.9 

76.0 

55.9 

55.2 

6 

FT 

87.7 

102.0 

75.9 

88.4 

7 

h+ 

118.0 

138.0 

103.0 

121.0 

8 

Ff 

161.0 

188.0 

142.0 

168.0 

9 

FT 

221.0 

259.0 

198.0 

234.0 

10 

H+ 

303.0 

355.0 

277.0 

327.0 

11 

fT 

417.0 

489.0 

387.0 

459.0 

12 

FT* 

574.0 

674.0 

543.0 

643.0 

13 

fT 

791.0 

929.0 

762.0 

903.0 

14 

H+ 

1091.0 

1281.0 

1071.0 

1269.0 

15 

FT 

1505.0 

2000.0 

1505.0 

2000.0 


source for the ion flux model, the Polar spacecraft, does provide information that allows for calculations of ion flux 
at high magnetic latitudes in the outer magnetosphere. The Polar spacecraft is in a polar orbit (2 Re x 9 Re) which 
crosses all magnetic latitudes. While the Polar spacecraft’s orbit does not directly sample the outer magnetosphere 
(> 9 Re), its data can be mapped along magnetic field lines to more completely fill the magnetosphere. 

Ion flux measurements used in creating the IFM database are from the Energetic Particles and Ion Composition 
(EPIC) Ion Composition Subsystem (ICS) on the Geotail spacecraft and the Comprehensive Energetic Particle and 
Pitch Angle Distribution (CEPPAD) Imaging Proton Sensor (IPS) from the Polar spacecraft. Table 1 provides a list 
of the EPIC/ICS ion energy channels. Energy channels for the CEPP AD/IPS instrument are listed in Table 2. 
Complete sets of data from 1995 through 2004 have been processed for the IFM database from the Geotail 
spacecraft. 1995 is when Geotail changed its orbit from sampling the distant magnetotail to its present orbit. The 
Polar ion flux measurements used date from the spacecraft’s launch in 1996 through 2004. 


2.1.1 Region Identification 

The first step in preparing the data for the IFM database generation process is to identify the phenomenological 
region (i.e., solar wind, magnetosheath, magnetosphere) each ion flux measurement is made in. A tool was created 
in IDL (Interactive Data Language) to allow for the graphical display of plasma moment data to identify boundary 
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Figure 5. Graphical interface to the geospace region identification tool. 



Figure 6. IFM Database Generation Process. 

crossing event times (Fig. 5). A record of each region boundary crossing event was made, allowing for the 
association of a flux measurement with a region. 


2.1.2 Database Generation 


There are three steps in creating the IFM database (Fig. 6). Step 1 is the creation of merged data files for each 
year that data is collected. These yearly files contain information about the magnetic activity, spacecraft position, 
ion flux measurements, and the phenomenological region (i.e., solar wind, magnetosheath, magnetosphere). Step 2 
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is the creation of the runtime database files for each region. It is during this process that the flux is mapped and the 
statistics are compiled for elements of space whose volumes are 1 Re x 1 Re x 1 Re. 

Generation of the model’s database is a computationally intensive process since it requires mapping all flux 
values in the Polar and Geotail data sets to values in the model database. Individual databases are required for each 
energy, Kp value, or other input parameter, requiring significant computer resources to generate the model. 
However, once the databases are complete, the runtime code is computationally efficient since the time consuming 
calculations take place during database generation. Database generation is implemented in three steps: 

Step 1: Create a database of “streamline” position points that trace out the drift path available to charged particles 
while conserving both the total energy and the first adiabatic invariant (the magnetic moment) as they propagate 
through the magnetosphere. The particle drift velocity V D is 


D 


= V 


VB 


+ V 


Curv 


+ V 


ExB 


(i) 


where 


Vyg = B-field gradient drift velocity 
v Curv = field curvature drift velocity 
V ExB = ExB drift velocity. 

The total drift velocity of a particle in and electric field E and magnetic field B is [Parks, 1991] 


where, 


ExB m 

v = + 

D ^2 


2 ^ 
v 1 
v„ + 


B 


qB 


V 


Bx- 


VB 


/ 


q = electric charge of the particle 

V || = V COS a = particle speed parallel to magnetic field line 
V_l = V sin# = particle speed perpendicular to magnetic field line. 


( 2 ) 


Tsyganenko B-Magnitude (nT) 
Zero Geodipole Tilt (+0.023 deg) 
Zgsm = 0 

Kp = 3 => Dst = -1 7.1 78 



Mcllwain Kp Dependent Geoetectric Potential 
Kp = 3 
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XMlN = -15 Re; PH!DIFF = 2kV 
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(a) (b) 

Figure 7. Magnetic and Electric Potential Models, (a) Tsyganenko magnetic field intensity |B| (nT) and (b) 
geoelectric potential (kV) in the ZGSM = 0 plane. These values will be used to compute an example set of 
streamlines shown in later figures. 6 
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The Tsyganenko geomagnetic field model [Tsyganenko, 1995, 1997; Tsyganenko and Stern 1996] is used in the 
current version of the model since it provides magnetic field values as a function of solar wind plasma parameters 
and geomagnetic disturbance values. Figure 7a shows Tsyganenko total field intensity |B| in the Zgsm = 0 plane. 
The electric field is 


E = -VO 


(3) 


ExB - Grad-Curve Drill Streamlines 
Tsyganenko + IGRF + Mcllwain 

1 50 keV Protons; 90 deg Pitch Angle; Kp = 3 
Backplane: Xstrt = -15 Re 

THET AMAX = 160 deg. PHITOL = 5 0; RDEL = 0 25;Kp = 3 
30. 
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Figure 8. Example Streamlines and Flux Observations, (a) Streamlines generated in Step 2 are projected onto 
the ZGSM = 0 plane and can be compared to (b) Geotail spacecraft flux data for all phenomenological regions 
(solar wind, magnetosheath, and magnetosphere) without streamline mapping. 


where O is the geoelectric potential in the equatorial plane, currently provided by the Mcllwain K p -dependent 
model [ Mcllwain , 1986]. This model includes both the convective and co-rotating contributions to the electric 
potential. Magnetic field lines are treated as equipotential lines, allowing the potential at the equatorial plane to be 
mapped to higher magnetic latitudes. Figure 7b shows an example output from the Mcllwain geoelectric potential 
model. The modeling approach used allows for the straightforward substitution of alternate geomagnetic field and 
geoelectric potential models if higher fidelity models are acquired or tests of particle loading into the database 
utilizing alternative models are desired. Higher fidelity global geoelectric field models will be evaluated for use and 
can be easily incorporated into the model. 


The total particle energy is 


Ejot — KEd KEf PE 


(4) 


where, 


KE d = kinetic energy due to particle drift 
(gradient, curvature, ExB) 

KE f = kinetic energy from particle’s spiral 

motion about the B field line 
PE = particle’s potential energy. 

Alternatively, the energy equation (4) can be expressed in terms of the magnetic moment and geoelectric potential as 
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( 5 ) 


"TOT 


= ^ m ( v D +v ii)+- uB+ q <I> 


with the magnetic moment given by 


1 2 

— mv I 

(6) 

The streamlines are created by initializing the start positions of an autonomous agent (AA) algorithm representing 
test particles. Streamlines are found by minimizing deviations in energy and magnetic moment in the AA algorithm 
along the particle drift trajectory. The AA actions are based upon the calculated environments for the magnetic and 
electric fields for a range of activity levels, conserving particle magnetic moment and energy at each step in the 
calculation. Figure 8a shows an example of streamlines created for the electric and magnetic. The density of 
streamlines used in the actual database generation process is much higher than shown in the figure, which used a 
smaller number of streamlines for clarity. The streamlines in this case compare favorably to the flux distributions in 
the data sets as shown in Figure 8b where EPIC/IPS proton flux in the 100-200 keV energy range (for 2 < K p < 4) 
is projected onto the equatorial plane. An example of the streamlines compared to satellite measurements of 
magnetospheric ion flux are shown in Fig. 9, demonstrating that the streamline contours generated analytically are in 
good agreement with the spacecraft flux measurement data (ignore the streamline artifacts in the data “hole”). 



Magnetosphere 

XY Flux Slice (protons/(cm A 2-sec-sr-MeV)] 
Average of all Z-values; 1 00 - 200 keV proto: 
(Kp 2*4; Includes Solar Event Particles) 
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= -17.15 


Xgsm (Re) 
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XY Flux Slice [protons/(cm A 2-sec-sr-MeV)] 
Average of all Z-values; 1 00 - 200 keV protons 
(Kp 2-4; Includes Solar Event Particles) 



(a) (b) 

Figure 9. Streamline Overlay on Magnetospheric Ion Flux Distributions, (a) Ion flux within the 
magnetosphere are projected onto the Zgsm = 0 plane, (b) Streamlines shown in Figure 2a are plotted over the ion 
flux distribution. 


Step 2: Create a database of pointers that allow for the rapid mapping of a satellite particle flux measurement to a 

streamline. This cross-referencing database divides the magnetosphere into spatial volume elements that are at a 
much finer resolution than used in the final, runtime database used to perform flux calculations. 


Step 3: Generate the runtime database, using spacecraft (Geotail, Polar) particle flux measurements. Use the region 
crossing database to associate each spacecraft measurement with a phenomenological region (solar wind, 
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(a) (b) 

Figure 10. Magnetospheric Ion Flux Derived from Streamline Mapping, (a) Ion flux data set without 
mapping projected onto the GSE Z=0 plane and (b) the same data used with the streamline mapping algorithm. 

magnetosheath, magnetosphere). Use the cross-referencing database to determine which streamline(s) to attach a 
given satellite measurement of magnetospheric flux to. Bin the data, based upon the magnetic activity index (e.g., 
K p ). Allow the flux to be mapped up and down the streamline for a relatively small distance, performing range- 
weighted averaging of the particle flux. 


From Roderer [1970] we get the expression for the omnidirectional flux at a point P, expressed as an integral 
over the equatorial flux: 


j 0 (cos a 0 )cos or 0 d(cos a 0 ) 


where 

a = particle pitch angle, 

J p = omnidirectional flux at the point P 
j 0 = directional equatorial flux 
B p = magnetic field at the point in question 
B 0 = magnetic field at the field line intersection 
with the equatorial plane 
B e = magnetic field at the field line intersection 
with the absorbing atmosphere. 

Use of equation (7) allows flux to be mapped along the geomagnetic field lines from the streamline point in the 
equatorial plane. Accumulation of flux values in volume elements from multiple streamlines produces a position 
dependent database that is used to compute the ion flux statistics within each volume element. Range weighting 
ensures that the average flux calculated for a given volume element will reflect those measurements closest to the 


1 — — 

B n v 


cos 2 O' 


bA 

J(B P ) = 4*/ \ 
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Electric Field (Vim) 


volume element. Figure 10 shows that the mapping takes a sparsely populated database and more fully populates 
the magnetosphere database, while maintaining the overall structure. 

Examples of output from the database generation process are shown in Figures 2 through 4. These flux maps 
show the proton flux as a function of the Kp magnetic activity index intervals. These are scatter plots of the data in 
the x-y plane, showing the average of all z- values. The data clearly orders up by Kp, especially in low energy plots 
(Fig. 2). These plots also demonstrate the sparseness of the data as Kp increases. As the energy increases, the flux 
decreases and the distinction between the solar wind, magnetosheath, and magnetosphere regions becomes less. As 
expected, the higher energy particles have free access to the magnetospheric regions modeled. 


2.1.3 High Energy Ion Analysis 
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Figure 11. Diagnostic output from high energy ion Lorentz force propagator. Starting conditions for this test 
case: Xgsm = -12 Re; Ygsm = +5Re; Zgsm = +lRe; Kp = 3.0; KE = 1 MeV; Pitch Angle = +20 deg. 

The streamline mapping algorithm described above has been successfully used in generating the databases for 
CRM. However, the guiding center approximation used for CRM will not work in cases of high kinetic energy. 
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where the particle’s gyroradius has the possibility of exceeding the magnetic field’s gradient scale length. It is for 
this reason that a test module was created which uses a full physics based approach that can handle non-adiabatic 
particle motion. The major components of the high energy ion module are: the Tsyganenko geomagnetic field 
model, the Mcllwain geoelectric potential model, and a Lorentz force particle mover. There are three different 
implementations of the particle mover’s integration scheme: Euler’s method, Runge-Kutta with adaptive step-size, 
and Bulirsch-Stoer with adaptive step-size. Test cases were exercised with all three integrators. Euler’s method is 
the simplest, but results indicated it lacked the accuracy required to propagate ions in the outer magnetosphere. 
Both the Runge-Kutta and Bulirsch-Stoer integrators performed well. 


Figure 1 1 shows diagnostic output from one of the test cases using the Runge-Kutta integrator. From the range 
versus elapsed time calculations, it is seen that the 1 MeV particle starting at the dusk side of the outer 
magnetosphere with a 20 degree pitch angle mirrors after about 5 seconds at a geocentric distance of about 5 Re. 
The particle then heads in the tailward direction where the motion becomes chaotic once it crosses into the neutral 
sheet region, where the gyroradius of the particle greatly increases. However, from this and other numerical 
experiments it was demonstrated that the gyroradius remained relatively small except for cases where the particle 
entered the neutral sheet region. This and other test results suggested that the streamline methodology described 
earlier can be used for the energy range of particles included in the IFM database (58.1 to 1560 keV). 


Sector Locations Flux Variation in Sector 2 



(a) (b) 

Figure 12. Flux Scaling. Sectors are identified in (a) that were used to obtain scaling laws for the flux (in units 
of protons/[cm2-sec-sr-MeV]) as a function of Kp. The flux variation as a function of Kp is given in (b) for 
Sector 2. Statistics of the flux variation for each Kp level is determined to provide a basis for estimating the 
probability of encountering extreme flux values and a method of scaling statistically significant scaling flux 
distributions at low Kp to the poor statistics in the high Kp cases (from Blackwell et al ., 2000a). 


2.1.4 Flux Scaling 

Flux measurements are sparse at large Kp (Kp > 6) over much of the magnetosphere. For this reason a “flux 
scaling” approach has been adapted from CRM [. Blackwell et al ., 2000]. Scaling laws are used to scale the average 
flux for a given Kp to the mean, 50%, and 95% levels for any Kp value. Figure 12(a) shows that a number of spatial 
sectors were chosen so that Kp scaling laws could be derived which cover the entire range of Kp values. Figure 
12(b) shows an example sector’s results. Unique scaling laws are derived for each sector. 
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2.1.5 Polar Spacecraft Data 


IFM’s output energy bins are tied to the Geotail energy bins (Table 1). The number and energy limits of data 
channels from the Polar spacecraft (Table 2) differ from Geotail’s, so the database generation software was modified 
so that ion flux from Polar is calculated for the Geotail energy channels. This is performed by performing a power- 
law fit to the spectrum for each Polar flux measurement (Fig. 13). The power-law representation of Polar data 
allows for the integration of the flux inside each of the Geotail energy bins. 

For each point in the Polar data the power law fit of the flux is calculated over all the energy bands. 

Flux = AE b , (8) 

where E is the energy, A = e a , and B = b. The least squares fitting coefficients are defined as 

b = nyYlnXi lnvd - Yflnxd Yflnyfl (9) 

nHlnxOMllnxi) 

a = YOnv) - bYOnx). (10) 

n 


Polar Average Flux vs. Energy 
Day 5.60971 of 1999 



Figure 13. Power-law fit to a Polar spacecraft flux measurement. 

where x is the energy and y is the flux. The flux for Polar in the Geotail energy band was determined from the 
power law fit. 


Flux= AtE2 ffl ±a -El ffl±1} l (11) 

(B+1)(E2 - El) 

where E2 and El are the upper and lower bounds of the energy band, respectively. A and B are defined above. 


III. Runtime Model 

The runtime code for CRM and IFM share a great deal in common in their overall structure (Fig. 14). However, 
significant modifications were made so that IFM can handle not just one, but multiple energy bands in calculating 
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Figure 14. Major software modules in IFM. 


flux. Figures 15 - 17 are examples of output from IFM. These are scatter plots of the proton flux predicted by IFM 
in the Zgsm = 0 plane for a range of Kp geomagnetic activity indices. Figure 18 shows the mission fluence 
calculated for each of the energy channels available in IFM, using the Chandra ephemeris as an example. These 
plots show the utility of mapping the data, since the code fills in large spatial gaps between flux values, while still 
maintaining the overall spatial distribution of the flux. The user can control the extent of the near-neighbors used in 
calculating flux at a point. Also available to the user are a variety of smoothing algorithms. 


IFM P2 CWpul 
Zgsm = 0 Kp = 1 





IFMPZOUpul 
Zgvn » 0 Kp • 3 



IFM P2 OutpU IFM P2 OJtfU 



Figure 15. Kp variation for P2 channel with spatial sector weighting. 


IFM PS Output IFM PS 0<p*l 



IFM P5 OJ$*l 
~0 Kp = 5 


m 


4 





CFM PS Output 

ZyvnxO. Kp = 7 ■— 



Figure 16. Kp variation for P5 channel with spatial sector weighting. 
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IFWP3 Output IFMPSOtfpUl iru F3 OJtfOt IFMPaO Upd 



Figure 17. Kp variation for P8 channel with spatial sector weighting. 


The user must decide whether to use the Kp-correlated database included with IFM vl.O or to use a Kp + ACE- 
correlated database (not yet available). 


III. Summary 

IFM is a radiation environment risk mitigation tool that provides magnetospheric ion flux values for varying 
geomagnetic disturbance levels in the geospace environment. IFM incorporates ion flux observations from the Polar 
and Geotail spacecraft in a single statistical flux model. IFM is an engineering environment model developed to 
predict the proton flux in the solar wind, magnetosheath, and magnetosphere phenomenological regions of geospace. 
The development of IFM has heavily leveraged experience gained in creating the ion flux model for the Chandra X- 
ray observatory.. This database allows for IFM output to be correlated with the geomagnetic activity level, as 
represented by the Kp index. 


■ P2_Fluence 

• P3_Fluence 
A P4_Fluence 

♦ P5_Fluence 

★ P6_Fluence 
P7_Fluence 

E P8_Fluence 
O P9_Fluence 

0 500 1000 1500 2000 2500 3000 

MET (days) 

Figure 18. IFM fluence (#/cm2-sr-MeV) for Chandra ephemeris. 
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Data Overview: X-Z Plane 
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Polar CEPPAD/IPS Energy Bands 
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Major Software Modules in IFM 
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Scaling laws are used to scale the mean, 50%, and 95% levels for 
any Kp value 
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Geotail’s 

- Power law fit made to the spectrum for each Polar measurement 

- Allows for integration of flux inside each of the Geotail energy bins 
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IFM Fluence Calculation 
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The Living With a Star Ion Flux Model (IFM) is a radiation environment risk 
mitigation tool 
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IFM is more flexible and will have a wider applicability than CR 


